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INTRODUCTION 


High  frequency  (HF)  skywave  field  strength  and  transmission  loss 
predictions  are  necessary  for  a  variety  of  reasons,  apart  from  any  theoretical 
interest  in  making  such  calculations.  For  example,  in  establishing  a  com¬ 
munication  circuit  we  must  have  a  preliminary  estimate  of  expected  median 
system  performance.  Also,  effective  management  of  frequency  assets  requires 
knowledge  of  the  effect  of  ionospheric  characteristic  variations  on  system 
performance.  These  variations  are  then  used  to  provide  updated  predictions  of 
optimal  operational  frequencies. 

For  real-time  frequency  management,  we  normally  begin  with  long-term 
predictions  of  expected  performance  based  on  median  values  of  input  param¬ 
eters,  i.e.,  fQF'2 ,  sunspot  number,  etc.  These  predictions  are  then  updated 
with  locally  measured  values  of  one  of  the  relevant  parameters  to  give  a 
real-time  prediction. 

In  circuit  planning  applications,  we  are  primarily  interested  in  expected 
hourly  median  behavior.  Sometimes,  these  median  predictions  are  augmented 
with  predictions  of  shorter  time-scale  behavior.  For  example,  in  some  field 
strength  predictions,  models  that  predict  fading  rates  and  amplitudes  due  to 
ionospheric  scattering  are  employed.  Th ’se  variations  about  a  median  value 
may  have  time  scales  on  the  order  of  minutes  or  even  seconds.  However,  these 
models  are  also  based  on  median  behavior  determined  from  extensive  statistical 
analysis  of  global  scattering  data.  In  general,  predictions  of  field 
strengths  on  less  than  an  hourly,  monthly  median  basis  require  sounding  of  the 
ionosphere  to  determine  its  current  state  and  to  monitor  its  short-term 
variations . 

This  report  deals  with  the  long-term  monthly  median  prediction  of  field 
strength  or  transmission  loss.  Consequently,  all  physical  parameters  will  be 
assumed  to  be  monthly  median  values.  Also,  we  will  be  concerned  exclusively 
with  earth-based  transmitters  and  receivers.  For  the  most  part  this  restricts 
the  frequency  range  to  the  HF  band  of  2-30  MHz.  Later,  we  will  discuss 
possible  propagation  modes  at  frequencies  greater  than  30  MHz. 
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FIELD  STRENGTH  PREDICTION 


When  developing  a  method  for  prediction  of  field  strength  or  transmission 
loss,  the  ultimate  objective  is  obviously  a  simple  and  accurate  method.  In 
many  instances,  the  method  should  also  be  fast  enough  for  use  in  real-time 
applications.  In  general,  there  are  two  approaches  that  can  be  followed  to 
produce  such  a  model.  One  is  to  fit  analytic  equations  to  the  experimentally 
determined  dependence  c'  transmission  loss  on  path,  time,  and  frequency.  The 
other  approach  is  to  estimate  the  total  transmission  loss  as  the  sum  of  a 
number  of  separate  mechanisms  for  energy  loss.  In  this  method,  each  term  is 
an  expression  deduced  either  from  theory  or  measurement  or  both.  Both 
techniques  have  advantages  and  disadvantages  that  must  be  recognized. 

The  former  method  is  usually  simpler  in  its  implementation  but  requires 
an  extremely  inclusive  data  base  to  insure  that  all  observable  trends  have 
been  included.  Even  in  instances  where  an  adequate  data  base  exists,  errors 
inherent  in  the  fitting  process  produce  corresponding  errors  in  predictions 
(Second  CCIR  Compucetr-based ,  1978). 

The  latter  method  i.  conceptually  more  elegant  and  enables  variations  to 
be  specified  in  a  physically  meaningful  manner.  However,  failure  to  include 
an  important  term  can  lead  to  significant  error  in  predictions.  Also,  the 
usual  tendency  is  to  produce  a  somewhat  more  complicated  model  that  may  not  be 
justified  in  terms  of  the  gain  in  accuracy  achieved  (Second  CCIR  Computer- 
based,  1978)  . 

The  two  field  strength  models  tested  and  described  in  this  report  include 
examples  of  both  of  the  above  methods.  The  first  model  we  describe  is  that 
used  in  the  program  HFBC84 ,  which  is  the  latest  version  of  LIL252-2.  This 
program  represents  a  simplified  treatment  of  the  methods  established  by  the 
International  Radio  Consultative  Committee  (CCIR)  in  its  original  interim 
report  (CCIR  Interim  Method,  1970),  and  in  later  updates  (Second  CCIR 
Computer-based ,  1978).  The  second  method  is  that  used  in  the  Naval  Ocean 
Systems  Center  (NOSC)  PROPHET  family  of  prediction  programs.  This  method  is 
based  on  an  empirical  model  developed  at  Forschungsinstitut  der  Deutschen 
Bundespost  (FTZ)  by  Beckmann  (1967)  and  further  refined  by  Damboldt  (1975). 

The  HFBC84  version  of  the  CCIR- recommended  field  strength  calculation 
method  is  used  here  because  it  compares  most  closely  with  the  philosophy  and 
intent  of  the  PROPHET  program  to  provide  relatively  fast,  simple  estimates  of 
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field  strength.  This  speed  is  achieved  in  PROPHET  by  using  the  FTZ  method  in 
conjunction  with  the  MINIMUF  algorithm  for  the  maximum  usable  fiequency  (MUF) 
determination.  HFBC84  achieves  its  speed  of  execution  through  its  con¬ 
siderable  simplification  of  the  field  strength  models  employed  as  compared  to 
its  parent  programs. 

A  necessary  input  to  any  HF  skywave  field  strength  prediction  scheme  is 
the  basic  MUF.  For  a  particular  time  and  circuit,  this  is  defined  as  the 
maximum  frequency  wave  that  can  propagate  by  ionospheric  refraction  alone.  We 
begin  by  describing  the  two  methods  for  predicting  this  parameter  used  in  the 
respective  field  strength  prediction  programs  mentioned  above.  The  first  is  a 
simplified  version  of  the  method  recommended  by  the  CCIR  and  used  in  HFBC84. 
Next,  we  outline  the  MINIMUF- 85  algorithm  used  by  PROPHET. 

MUF  PREDICTIONS 
Simplified  CCIR  Method 

For  applications  where  speed  is  not  essential  and  adequate  computer 
memory  is  available,  a  method  for  predicting  the  median  F2  layer  basic  MUF 
based  on  numerical  maps  of  the  ionospheric  characteristics  f|^F2  and  M(3000)F2 
has  been  developed  (Report  to  the  Second  Session,  1984).  Here  ^^^2  is  the  F2 
layer  ordinary  wave  critical  frequency.  It  represents  the  maximum  frequency 
wave  reflected  back  to  the  earth  for  vertical  incidence  on  the  ionosphere. 
Numerically  it  is  proportional  to  the  square  root  of  the  peak  electron  density 
of  the  layer.  M(3000)F2  is  a  path  length  factor  that  represents  the 
multiplier  for  the  fgF2.  The  product  of  f^F2  and  M(3000)F2  gives  F2(3000)MUF, 
the  F2  layer  MUF  for  a  3000-km  path.  Both  of  these  characteristics  are  scaled 
from  vertical  incidence  ionograms . 

For  great  circle  path  lengths,  D(km) ,  less  than  4000  km,  the  F2- layer  MUF 
is  obtained  by  first  determining  f^FZ  and  M(3000)  at  the  path  midpoint  for  12- 
month  running  mean  sunspot  numbers  of  R^^^  =  0  and  R^^  =  100.  The  values  of 
these  characteristics  corresponding  to  the  required  R^^  index  are  found  by 
interpolation  if  the  required  value  is  between  0  and  100  and  by  extrapolation 
if  the  value  is  greater  than  100.  If  the  required  Rj^2  value  is  greater  than 
150,  the  value  at  150  is  used. 
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We  then  have 


F2(0)MUF  =  f  F2  +  f„/2 
o  H 

F2(4000)MUF  =1.1  •  f  F2  •  M(3000)F2 

o 

where  f  is  the  electronic  gyrofrequency  at  the  midpoint  of  the  path.  The 
F2- layer  MUF  for  the  path  distance  D  is  then  found  by  interpolation  according 
to  (Report  to  the  Second  Session,  1984), 

F2(D)MUF  =  F2(0)MUF  +  [F2(4000)MUF  -  F2<0)MUF]  •  M(D) 

where 

M(D)  =  1.64  •  10'^  •  ,  0  <  D  <  800 

or 

-14  4  -10  3  -72  -4 

M(D)  =  1.26  •  10  .  D  -  1.3  .  10  •  D  +  4.1  •  10  •  D  -  1 . 2  •  10  •  D , 

800  >  D  <  4000  . 

For  great  circle  path  lengths  greater  than  4000  km,  a  simple  control 
point  method  is  employed.  Control  points  are  taken  at  points  2000  km  from 

each  end  of  the  path,  and  at  these  points  the  value  of  F2(4000)MUF  is 

determined,  interpolating  in  required.  The  lower  of  these  two  values 

is  taken  as  the  median  F2- layer  MUF  for  the  path. 

For  all  path  lengths,  the  F2- layer  MUF  is  calculated  according  to  the 
above  method.  For  paths  between  0  and  4000  km,  the  median  basic  MUF  of  an 
E-layer  mode  is  also  determined.  This  is  done  by  first  determining  the 
E- layer  critical  frequency,  at  the  path  midpoint  for  a  path  of  2000  km  or 

less,  or  at  points  1000  km  from  each  end  for  a  path  between  2000  and  4000  km. 
ac  these  joints  the  solar  zenith  angle  x  is  determined  and  the  f^E  is  given  by 
(Report  to  the  Second  Session,  1984) 

f^E  =  0.9  •  [(180  +  1.44  .  *  cos  MHz 

where 

X'  =  X.  0  >  X  >  80° 


80°  <  X  <  116° 


X'  =  89.907  ,  X  >  116” 
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The  E- layer  basic  MUF  for  the  path,  E(D)MUF,  is  then  given  by 

E(D)MUF  =  f  E  •  sec  i,,.  (1) 

o  iiU 

where  f^E  is  the  midpoint  critical  frequency  for  paths  less  than  2000  km  and 
the  lesser  of  the  two  values  for  paths  greater  than  2000  km.  Here  is  the 

angle  of  incidence  of  the  ray  at  a  height  of  110  km  assuming  mirror  reflection 
in  the  ionosphere. 

For  most  applications,  the  F2  layer  controls  the  propagation  charac¬ 
teristics  on  a  given  path,  and  the  former  method  based  on  numerical  maps  is 
required.  The  number  of  numerical  maps  required  to  specify  the  value  of  f^F2 
and  M(3000)F2  for  each  hour  of  every  month  at  =  0  on  R^2  “  iO®  is 
extremely  large.  Add  to  this  the  fact  that  the  median  behavior  of  these 
characteristics  can  become  quite  complicated,  and  we  see  that  the  total  number 
of  coefficients  required  for  this  method  makes  it  impractical  for  micro¬ 
computer  application.  Consequently,  techniques  specifically  developed  for  use 
on  microcomputers  are  now  used  in  many  applications.  Next,  we  describe  one  of 
the  faster  and  shorter  methods,  MINIMUF-85. 

MINIMUF-85 

MINIHUF-85  is  the  latest  version  of  the  MINIMUF  algorithm,  which  has  been 
described  many  times  in  other  documents  (Rose  and  Martin,  1978;  Sailors  et 
al . ,  1986).  All  PROPHET  systems.  Advanced  Classic  PROPHET,  APES,  TDA,  etc., 
use  the  MINIMUF  algorithm  to  predict  monthly  median  values  of  the  F2- layer 
MUF.  As  mentioned  earlier,  the  advantage  of  the  MINIMUF  method  is  its 
comparably  fast  execution  time,  its  simplicity,  and  its  ease  of  transport, 
which  makes  it  particularly  suitable  for  real-time  and  field  use.  Sailors  et 
al .  (1981)  found  it  to  be  extremely  accurate  considering  the  simplicity  of  the 
method. 

in  MINIMUF-85,  the  F2- layer  MUF  for  a  path  of  length  D(km)  is  given  by 

F2(D)MUF  =  f  F2  •  M(D) 
o 

where 

f^F2  =  (6  +  Al  •  Jcos  xeff)^^^  (2) 
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and 

M(D)  =  A2  •  A3  •  A4  .  {1  +  2.5  •  sin(2.5  •  . 

Here,  Al  and  A2  are  linear  functions  of  monthly  median  sunspot  number,  A3  is  a 
sixth-order  Fourier  series  based  on  season,  and  A4  is  a  function  of  time,  the 
form  of  which  differs  between  night  and  day.  The  great  circle  path  length 
corresponding  to  D,  in  radians,  is  4'.  The  model  uses  a  control  point  method 
for  paths  greater  than  4000  km  similar  to  that  described  earlier. 

In  equation  2,  i-S  an  effective  solar  zenith  angle  and  cos  is 

modeled  as  the  "response  of  a  linear  first-order  system  driven  by  the  actual 
cos  x"  follows  (Rose  and  Martin,  1978), 

^  (cos  xsff)  +  cos  xcff  =  cos  X  • 

An  effective  zenith  angle  is  used  to  simplify  the  modeling  technique  since  the 
F2-layer  variation,  unlike  the  E  layer,  is  not  well  represented  as  a  function 
of  zenith  angle  alone.  By  using  an  effective  zenith  angle  and  fitting 
constants  in  the  model  to  actual  measured,  oblique-path  MUFs,  the  effect  of 
the  other  dependencies  of  the  F2  layer  can  be  accounted  for  without  explicit 
inclusion  in  the  model. 

As  of  this  writing,  MINIMUF-85  does  not  include  a  determination  of  the 
E- layer  MUF.  We  anticipate  in  the  future  a  calculation  similar  to  that  used 
in  the  CCIR  program  will  be  included  for  short  (<2000  km)  paths. 

A  statistical  comparison  of  the  accuracy  of  the  above  methods  for 
predicting  MUFs  is  shown  in  table  1.  These  results  are  based  on  a  comparison 
of  predicted  values  with  over  13,000  measured  oblique  sounder  MUFs.  Table  1 
shows  that  the  two  methods  of  MUF  predictions  are  equally  accurate  for  the 
most  part  and  can  be  expected  to  predict,  on  the  average,  approximately  equal 
MUFs  as  input  to  the  individual  field  strength  prediction  programs.  A  more 
complete  comparison  of  the  two  MUF  prediction  methods  is  presented  in  a  NOSC 
technical  report  (Roy  and  Sailors,  1987). 
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Table  1.  Summary  of  comparison  of  MINIMUF-85  and 
HFBC84  MUF  predictions. 


Conditions 

MUF85 

HFBC84 

Population  size 

13,054 

13,054 

Average  residual 

1.282 

1.165 

RMS  residual 

4.579 

4.665 

MAE  residual 

3.451 

3.472 

Average  relative  residual 

.051 

.059 

RMS  relative  residual 

.239 

.242 

MAE  relative  residual 

.180 

.179 

Average  absolute 
relative  residual 

.186 

.196 

Standard  error  of  estimate 

3.971 

3.885 

Correlation  coefficient 

.819 

.827 

Next  we  present  a  description  of  the  two  field  strength  models.  We  begin 
with  the  methods  used  in  HFBC84  and  then  present  the  PROPHET  field  strength 
model.  A  relatively  detailed  description  of  the  models  will  be  given  since 
these  serve  to  highlight  the  various  mechanisms  through  which  energy  is 
reduced  (or  increased)  for  a  wave  propagating  in  the  ionosphere. 

FIELD  STRENGTH  CALCULATIONS 

HFBC84  Method 

Path  Lengths  <7000  km.  The  method  used  for  these  path  lengths  assumes 
great  circle  geometrical  ray  propagation  with  reflection  from  E  or  F2  layers 
in  a  horizontally  stratified  ionosphere.  The  restriction  to  path  lengths  less 
than  7000  km  is  due  to  the  recognition  that,  while  the  ionosphere  is 
dominantly  horizontally  stratified  due  to  the  nature  of  the  ionization 
process,  there  is  also  a  high  degree  of  variability  and  fluctuation  in  the 
ionization  levels  present  at  all  times.  For  a  wave  passing  through  this 
variable  medium,  the  majority  of  the  energy  is  refracted  into  the  great  circle 
path.  However,  scattering  due  to  variability  in  the  form  of  layer  tilts  or 
electron  density  irregularities  causes  some  propagation  out  of  this  path.  For 
long  path  lengths,  the  wave  may  pass  into  the  ionosphere  four  or  more  times 
with  correspondingly  more  of  its  energy  scattered  out  of  the  great  circle 
path,  thus  decreasing  the  received  field  strength.  Likewise,  for  a  given 
great  circle  transmission  path,  energy  from  adjoint  rays  that  would  not 
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nor^nally  be  expected  to  contribute  may  be  scattered  to  the  receiver,  causing 
an  increase  in  the  received  field  strength.  As  currently  modeled,  the  ray 
path  geometry  method  of  calculation  is  unable  to  account  for  this  non- great 
circle  propagation,  and  other  methods  are  used  for  longer  path  lengths  where 
these  effects  become  increasingly  important. 

The  first  step  in  calculating  the  field  strength  is  determining  the  MUF 
for  each  mode  to  be  included.  The  minimum  hop  F2-  and  E- layer  MUFs  are 
calculated  as  outlined  above.  The  MUF  for  higher  order  modes  is  found,  in  the 
case  of  the  E  layer,  by  determining  the  angle  of  incidence  for  a  given  mode  at 
a  height  of  110  km.  The  MUF  is  then  found  by  multiplying  the  secant  of  this 
angle  by  the  midpath  f^E  for  paths  less  than  2000  km  or  by  the  minimum  f^E  at 
the  two  control  points  for  paths  greater  than  2000  km. 

For  F2- layer  modes,  the  MUF  is  found  by  determining  the  angle  of 
incidence,  assuming  mirror  reflection,  at  a  height  given  by  the  Shimazaki 
equation  (Second  CCIR  Computer-based ,  1978) 


_  1490 

M(3000)F2 


-  176  (km) 


where  M(3000)  is  calculated  at  the  position  of  the  minimum  fQF2  at  the  two 

control  points  for  paths  greater  than  4000  km  and  at  the  midpath  for  paths 

less  than  4000  km.  The  secant  of  this  angle  is  then  multiplied  by  the 

minimum  f  F2 ,  or  midpath  f  F2 ,  to  determine  the  mode  MUF.  Five  modes  are 
o  o 

calculated  for  both  the  E  and  F2  layers. 

Once  the  MUF  for  a  mode  is  determined,  the  median  field  strength  (in  dB 
above  IpV/m)  for  that  mode  is  given  by 

Et3  =  136.6  +  20  log  f  -  -  L.  -  -  1^  -  -  7.3  .  (3) 

Here  136.6  -i-  20  log  f  represents  the  reference  field  at  1  km  for  a  wave  of 
frequency  f  (MHz)  for  1  kW  input  to  an  isotropic  antenna  in  free  space.  is 

the  actual  power  of  the  transmitter  in  dB  relative  to  1  kW  and  is  the  gain 
of  the  transmitting  antenna  in  dBi  in  the  direction  of  the  ray  for  the  mode. 
The  various  loss  terms  in  equation  3  will  be  dealt  with  separately. 

In  general,  for  paths  less  than  7000  km,  the  field  strength  values  of  the 
two  strongest  F2- layer  modes  and  the  strongest  E- layer  mode  are  determined 
according  to  equation  3.  The  resultant  field  strength  from  these  modes  is 
obtained,  assuming  uncorrelated  phases,  by  power  addition.  Since  multiple- 
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hop,  E- layer  modes  suffer  substantial  deviative  absorption,  these  modes  are 

not  considered  beyond  4000  km  {Report  to  the  Second  Session,  1984).  For  these 

distances,  only  multiple-hop,  F2-layer  modes  are  considered. 

a.  Basic  Free  Space  Loss,  The  major  source  of  energy  reduction  of 

the  wave  is  usually  the  free  space  loss  resulting  from  the  geometrical 

spreading  of  energy  as  the  radio  wave  progresses  away  from  the  transmitter. 

As  a  first  approximation,  for  a  one -hop  path,  the  CCIR  assumes  that  the  earth 

and  ionosphere  are  both  flat  and  that  reflection  is  specular  {CCIR  Interim 

Method,  1970) .  In  this  approximation,  the  energy  density  diminishes  as  the 

inverse  square  of  the  ray  path  distance  (no  focus  gain  at  these  path  lengths) . 

Now  for  an  isotropic  antenna  in  free  space  transmitting  P  watts  of  power, 

2 

the  power  flux  density  at  a  distance  D  (km)  is  P/4>rD  .  For  an  isotropic 

2 

receiving  antenna  in  free  space,  the  effective  receiving  area  is  X  /4?r,  where 
A  (km)  is  the  free -space  wavelength  of  the  incident  wave.  Thus,  the  total 
power  received  by  an  isotropic  antenna  excited  by  an  isotropic  antenna  at  a 
distance  D  is 


The  basic  free -space  transmission  loss  in  this  case  is  given  by  the  ratio  of 
power  transmitted  to  power  received,  or  in  dB, 


Sf  =  10  log 


Pf— 1 

t47rDj 


=  20  log  ^ 


This  can  be  written  as 


=  21.98  +  20  log  D  -  20  log  A  . 

Writing  A  in  terms  of  the  frequency  f  (MHz)  this  becomes 

=  32.44  +  20  log  D  +  D  log  f  .  (4) 

This  loss  is  due  to  spatial  expansion  of  the  wave  as  it  leaves  the  source  and 
it  is  present  in  any  application.  In  the  case  of  an  earth-based  transmitter, 
the  range  is  the  slant  path  length  of  the  ray  through  the  ionosphere.  For  an 
equivalent  triangular  path,  this  slant  path  range  is  given  by 
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d  =  2R 


cos [a  +  y 


for  a  1-hop  path.  For  an  arbitrary  n-hop  path  this  becomes 

D 


n 


sin" 


d  =  2R  I 


2R 


r, 

,  D  ] 

cos  1/ 

^  2rJ 

where  D  is  the  great  circle  hop  length  for  an  n-hop  mode  and  R  is  the  earth's 
radius.  Here  A  is  the  vertical  radiation  angle  for  the  mode  and  is  given  by 


A  =  tan 


-1 


R 


R  +  h 


,  cosec 


D_‘ 

2rJ  ■ 


where  h'  is  the  equivalent  mirror  reflection  height  in  the  ionosphere 
determined  for  the  mode.  For  E- layer  modes  h'  is  taken  to  be  110  km.  For  F2- 
layer  modes  h'  is  determined  at  the  control  point  with  the  minimum  value  of 
f^F2,  or  the  midpath  point  for  paths  less  than  4000  km.  It  is  given  by 


358 


(11  -  100a)  18.8  - 


320 

vS 


]  +  aD[o.03  + 


km 


or  500  km,  whichever  is  less.  Here  a  is  given  by 


1 

M(3000)F2 


0.24 


or  0.04,  whichever  is  greater.  X  is  the  ratio  of  F2- layer  critical  frequency 
to  E- layer  critical  frequency,  or  2,  whichever  is  larger. 

For  a  fixed  distance  D,  equation  4  shows  that  the  losses  increase  with 
frequency.  This  is  an  artifact  of  the  assumption  of  the  receiving 
characteristics  of  an  antenna  which  is  frequency  dependent.  However,  the 
free -space  loss  is  dependent  on  distance  only,  and  this  term  will  cancel  with 
a  similar  term  in  the  reference  field  and  will  not  appear  in  the  final  field 
strength  expression. 

b.  Ionospheric  Absorption  Loss,  L^.  The  next  most  important  field 
strength  reduction  mechanism  is  absorption  losses  suffered  by  the  wave  as  it 
propagates  through  the  ionosphere.  These  losses  come  about  through 
collisions,  mainly  between  electrons  and  neutral  particles  at  D- layer  heights 
and  electrons  and  positive  ions  at  E-  and  F-layer  heights. 
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As  an  electromagnetic  wave  travels  through  a  plasma  such  as  the 
ionosphere,  the  electrons  in  the  plasma  are  accelerated  into  motion  by  the 
electric  field  vector  of  the  wave.  (Here  we  are  assuming  that  at  HF  radio 
frequencies  the  motion  of  the  much  more  massive  ionic  component  of  the  plasma 
can  be  ignored  relative  to  the  electrons.)  In  the  absence  of  collisions  with 
other  particles  of  the  plasma,  the  accelerated  electrons  would  reradiate  their 
acquired  energy,  which  would  be  restored  to  the  passing  wave.  However, 
because  the  electrons  have  a  finite  mass,  this  radiated  energy  is  delayed  in 
time  with  respect  to  the  exciting  wave.  This  causes  a  phase  difference  to 
develop  between  the  two  energy  components,  resulting  in  a  change  in  the  phase 
velocity  of  the  total  wave  and,  hence,  in  the  index  of  refraction.  If  we 
include  the  effects  of  collisions  on  the  electrons,  the  above  picture  is 
significantly  altered.  The  electrons,  set  into  motion  by  the  wave,  collide 
with  other  particles  and  lose  some  of  their  acquired  kinetic  energy  to  these 
particles.  This  results  in  an  overall  increase  in  the  effective  temperature 
of  the  plasma.  The  energy  that  goes  into  heating  the  plasma  is  not  available 
to  be  returned  to  the  passing  wave  as  described  above  in  the  "no  collision" 
case.  In  this  way  energy  is  "absorbed"  from  the  wave  by  the  plasma  (Davies, 
1969)  . 

Mathematically,  the  absorption  process  in  the  ionosphere  is  described  by 
the  imaginary  part  of  the  complex  refractive  index,  x-  terms  of  x  we 

define  the  absorption  coefficient,  K,  as  the  imaginary  part  of  the  complex 
wave  vector. 


where  w  =  2nf  and  c  is  the  speed  of  light  in  vacuum. 

If  we  initially  ignore  the  earth's  magnetic  field  and  assume  an  isotropic 
electronic  plasma,  the  dispersion  relation  is  given  by  Booker  (1984) 


2 

2  2 
n^  =  (M  -  ix)  =  1  -  - 

U) 


where  ~  is  the  electronic  plasma  frequency.  Then,  the  absorption 

o 

coefficient  can  be  written  as 
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K  = 


Ni^ 


2e  ra/ic  2  2 

o  w  +  1/ 


(5) 


where  e  is  the  charge  of  the  electron;  is  the  free-space  permittivity,  m  is 
the  electronic  mass,  N  is  the  electron  density,  and  i/  is  the  effective 
collision  frequency. 

For  an  HF  wave  of  frequency  f^  launched  vertically  into  an  isotropic 
ionosphere,  the  total  path  absorption  is  given  by 


L(f  ) 

V 


2  /  Kdh 

h 


or 


L(f^) 


; 

h 


Nfh'")  i/(h')  dh' 

^  [IP]! 


where  is  a  constant,  h^  is  the  height  of  the  bottom  of  the  ionosphere,  and 
h^  is  the  height  of  reflection. 

Equation  5  shows  that,  for  a  fixed  operating  frequency,  the  absorption 
coefficient  can  become  large  if  the  product  Ni/  becomes  large  or  if  the  index 
of  refraction  becomes  small  while  Ni/  remains  finite.  These  two  limits  of 
equation  5  are  used  to  artificially  separate  the  absorption  suffered  by  an 
electromagnetic  wave  into  two  parts,  nondeviative  and  deviative. 

At  D- layer  heights  the  local  index  of  refraction  is  approximately  unity 
at  HF  frequencies  and  the  wave  suffers  little  refraction.  However,  at  these 
heights  the  product  N:/  can  become  large,  resulting  in  a  maximum  of  the 
absorption  coefficient.  This  is  referred  to  as  nondeviative  absorption.  In 
this  case  the  absorption  coefficient  is  given  by 


K  = 


8n  me  c 
o 


Ni/ 


f  ^  + 

V 


(6) 


If  the  wave  frequency  f^  is  less  than  the  maximum  critical  frequency  (E 
or  F  layer)  of  the  ionosphere,  the  wave  will  propagate  into  the  ionosphere 
until  it  reaches  a  height  at  which  the  plasma  frequency  is  equal  to  the  wave 
frequency.  At  this  level  n  =  o  and  the  wave  is  reflected.  For  usual 
operative  HF  frequencies,  this  occurs  in  the  E-  or  F- layer  region  of  the 
ionosphere . 
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Equation  5  shows  that  at  the  level  fi  =  o  the  so-called  deviative 
absorption  becomes  infinite  (since  Nt/  >  o  at  all  ionospheric  heights) .  This 
nonphysical  result  is  due  to  the  fact  that  the  geometrical  optics  approxima¬ 
tion  upon  which  the  above  argument  is  based  is  not  valid  near  the  turning 
point  for  the  wave. 

At  heights  near  the  reflection  level,  the  wave  suffers  considerable  group 
retardation.  Consequently,  the  wave  spends  a  long  time  in  the  absorbing 
region,  exciting  an  increased  number  of  collisions  and  causing  the  total 
absorption  to  increase  despite  the  low  neutral  particle  density.  At  these 
heights,  the  collisions  involve  positive  ions  instead  of  neutral  particles  and 
the  relations  become  more  complex,  so  we  will  not  pursue  the  topic  here. 

In  summary,  a  vertical  wave  propagating  through  the  ionosphere  suffers 
energy  loss  through  the  absorption  of  the  energy  by  collisions  with  neutral 
particles  at  D- layer  heights  (nondeviative  absorption)  and  by  collision  with 
positive  ions  (deviative  absorption)  at  E-  and  F-layer  heights.  Inspection  of 
equation  5  also  shows  that  this  absorption  is  approximately  inversely 
proportional  to  the  square  of  the  operating  frequency  for  nondeviative 
absorption. 

So  far  we  have  ignored  the  effect  of  the  earth's  magnetic  field.  Since 
this  field  greatly  affects  electromagnetic  wave  propagation  in  the  ionosphere, 
especially  as  the  wave  frequency  approaches  the  gyrofrequency ,  we  are  not 
surprised  to  encounter  its  effect  in  ionospheric  absorption.  In  fact,  a  more 
complete  analysis  shows  we  can  account  for  the  magnetic  field  in  our  earlier 
results  by  including  the  longitudinal  component  of  the  gyrofrequency  in  our 
initial  absorption  expression,  equation  5,  in  the  following  way 


K  =  a 


(7) 


where  a  is  a  constant  and  fj^  is  the  longitudinal  component  of  the  electronic 
gyrofrequency.  We  can  include  the  earth's  magnetic  field  to  this  level  of 
approximation  (quasi-longitudinal  approximation)  in  all  of  our  analysis  in  the 
same  way. 

Also,  in  our  analysis  above,  we  have  been  interested  in  vertical 
propagation.  Clearly,  we  wish  to  generalize  to  oblique  incidence 
applications.  This  is  accomplished  by  using  the  Martyn  absorption  theorem 
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(Davies,  1969).  This  theorem  states  that  for  a  plane-parallel,  horizontally 
stratified  ionosphere  in  the  absence  of  the  earth's  magnetic  field,  the 
absorption  L  experienced  by  a  wave  of  frequency  f^j^  incident  at  an 

oblique  angle  i^  on  the  base  of  the  ionosphere  is  related  to  that  experienced 
by  a  wave  at  vertical  incidence  of  frequency  f^  by 

L(f  )  =  L(f  )  cos  i  .  (8) 

Ob  V  o 

Here  the  waves  of  frequency  f^^  and  f^  are  equivalent  in  the  sense  that  their 
true  heights  of  reflection  are  equal.  Under  these  conditions,  the  secant  law 
states  that  (Davies,  1969) 


f  =  f  sec  i 
ob  V  o 


Substitution  of  this  into  equation  8  shows  that 


i 


L(f  )f  ^ 

^^^ob^  =  ~7~f~  ^ 

ob 

This  equation  is  exact  for  a  plane-stratified  ionosphere  with  no  magnetic 
field.  When  the  curvature  of  the  earth  and  the  magnetic  field  are  introduced, 
the  transformation  is  no  longer  exact.  However,  the  quasi -longitudinal 
approximation  remains  valid.  So  for  a  curved  earth  with  magnetic  field,  we 
can  approximate  the  nondeviative  absorption  loss  experienced  by  the  oblique 
wave  as 


(fob  "  fj' 


sec  1 


(10) 


If  we  further  assume  that,  for  the  most  part,  in  the  HF  band  we  have  f^  »  i/ , 
we  can  use  equation  7  to  give 


L(f 


Ob' 


A.  sec  i 
1  o 


(11) 


where  is  a  measured  index  of  the  absorption  experienced  by  the  vertical 
wave  of  frequency  f^,  which  is  equivalent  to  f^^  in  the  above  sense. 

For  the  HFBC84  program,  a  semiempirical  approximation  to  the  above  result 
has  been  developed.  Guided  by  the  form  of  equation  11,  measurements  on  a 
series  of  oblique  paths  have  led  to  the  following  approximation  for  the 
nondeviative  absorption  loss  term 
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L. 

1 


where 


677.2  sec 


i  n 

I  I- 

10.2  j=l  ^ 


n  =  number  of  hops  for  the  mode 


(12) 


i  =  angle  of  incidence  at  100  km 

f  =  operating  frequency  (MHz) 

f  =  full  gyrofrequency  at  100  km 
H 

and 

1  3 

=  (1  +  .0037  •  R^2)(‘^os  ®®1  x)j  where  R^2  1®  the  12-month  running 

mean  sunspot  number  and  Xj  Is  the  solar  zenith  angle  at  the  point  of  passage 

through  a  height  of  100  km  for  each  hop. 

At  night,  the  ionization  density  of  the  layers  decreases  and  the 

absorption  does  not  cease  but  falls  off  to  a  small  residual  value.  To 

approximate  this  behavior,  the  value  of  I .  is  restricted  to  be  greater  than 

0.1.  This  assumption  agrees  with  the  work  of  Wakai  (1961)  on  nighttime 

absorption  (Lucas  and  Haydon,  1966). 

Finally,  we  note  that  the  absorption  expression  of  equation  12  is  meant 

to  describe  the  median  level  that  can  be  expected  to  exist.  It  is  not  meant 

to  be  used  during  times  of  anomalous  solar  activity,  i.e.,  solar  flares,  when 

the  high  X-ray  flux  levels  emitted  by  the  sun  can  cause  extreme  levels  of  D- 

region  absorption.  As  we  mentioned  earlier,  since  nondeviative  absorption  is 

inversely  proportional  to  the  square  of  the  operating  frequency,  this 

increased  level  of  absorption  can  cause  the  minimum  available  operating 

frequencies  to  rise  to  the  level  where  it  actually  reaches  the  MUF.  During 

this  "shortwave  fade"  period,  no  skywave  is  available  for  communication. 

Behavior  of  this  sort  must  be  included  independently  in  separate  absorption 

models  derived  to  handle  these  effects. 

c.  Ground  Reflection  Loss,  L  .  For  multiple-hop  paths,  the  third  major 

& 

source  of  energy  loss  is  that  due  to  intermediate  ground  reflection.  Assuming 
the  wave  incident  at  the  ground  is  randomly  polarized  with  energy  equally 
distributed  in  horizontal  and  vertical  polarization,  the  loss  is  given  by 
{CCIR  Interim  Method,  1970) 


15 


L 

g 


10  log 


(13) 


where  and  are  the  vertical  and  horizontal  polarization  Fresnel 
reflection  coefficients,  respectively.  They  are  defined  as  the  ratio  of  the 
magnitude  of  the  electric  vector  in  the  reflected  wave  to  that  of  the  incident 
wave.  R^  and  R^  are,  in  general,  complex  and  are  given  by  the  Fresnel 
formulae 

2  22  1/2 

R^  =  n  sin  -  (n  -  cos  /3)  ' 

2  2  1/2 

P  _  sin  B  -  (n  -  cos  B)  ^ 

^  ^  /  2  2  ..1/2 

sin  ^  +  (n  -  cos  p) 

where  is  the  angle  of  the  reflected  wave  at  the  earth  and  n  is  the  complex 
refractive  index  of  the  earth. 

In  the  HFBC84  program,  these  calculations  are  not  performed  since  it  has 
been  determined  that,  to  an  acceptable  level  of  accuracy,  a  value  of  2  dB  can 
be  assigned  for  each  ground  reflection.  Thus,  we  have 


L 

g 


0,  for  I'hop  mode 
2  dB,  for  2 -hop  mode 
4  dB,  for  3 -hop  mode 
^  (n  -  1)  •  2,  for  n-hop  mode. 


(14) 


d.  Auroral  Absorption,  1^.  The  absorption  processes  in  auroral  regions 
are  somewhat  different  than  those  encountered  at  lower  latitude.  At  high 
latitude  (north  and  south) ,  the  principal  ionization  process  is  particle 
precipitation.  The  particle  flux  from  the  sun  is  carried  along  the  earth's 
magnetic  field  lines  and  deposited  at  ionospheric  heights  in  extremely  complex 
patterns  at  the  earth's  magnetic  poles.  This  flux  of  particles  causes 
ionization  through  both  collision  and  radiation  processes.  While  there  may  be 
some  solar  control  of  the  ionization  level  at  these  latitudes  (i.e.,  through 
X,  the  solar  zenith  angle  as  in  equation  12)  particle  precipitation  is  the 
major  influence  in  polar  regions. 

The  method  for  inclusion  of  these  processes  in  the  CCIR  program  is 
through  a  table  of  experimentally  determined  values  which  are  added  to  the 
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resultant  field  strength  (Report  to  the  Second  Session,  1984).  This  table 
contains  dependence  on  path  length,  time,  season,  and  geomagnetic  latitude. 

The  table  also  includes  some  small  low-  and  mid- latitude  corrections  used  to 
include  effects  not  explicitly  included  in  the  loss  terms  of  equation  3.  In 
general,  the  values  in  the  table  peak  at  polar  latitudes  with  values  from  3  to 
10  dB. 

The  next  loss  mechanism  given  in  equation  3  is  the  over- the -MUF  loss,  L  . 

m 

We  postpone  discussion  of  this  term  until  we  have  introduced  the  long  path 
model  used  in  HFBC84. 

The  numerical  factor  of  7 . 3  subtracted  from  equation  3  is  included  to 
account  for  "those  effects  of  skywave  propagation  not  otherwise  included  in 
this  fast,  simple  method"  (Report  to  the  Second  Session,  1984).  It  is  subject 
to  change  as  further  comparison  to  measured  data  is  completed. 

The  calculations  outlined  above  are  repeated  for  each  successive  higher 
order  F2  mode  until  a  peak  in  the  field  strength  is  found.  Then,  the  two 
highest  field  strength  modes  are  selected. 

For  E-region  modes,  the  lowest  order  mode  is  a  IE  for  ranges  of  0-2000  km 
and  2E  for  2000-4000  km.  For  paths  less  than  4000  km,  the  highest  field 
strength  mode  is  found  according  to  the  above  methods  and  the  result  is  added, 
as  power,  to  the  two  strongest  F2  layer  modes.  The  result  is  the  predicted 
field  strength  for  that  frequency  and  hour. 

Path  Lengths  Between  7000  and  9000  km.  For  paths  between  7000  and  9000 

km,  the  ray  path  geometry  procedure  for  F2-layer  modes  is  used  in  conjunction 

with  the  method  for  long  paths  (>9000  km)  to  be  described  below.  The  results 

from  the  two  methods,  E  and  E  are  combined  as  follows 

ts  t  r 


E  . 
ti 


D_-_7^0^ 

ts  2000  ^  t^ 


(15) 


where  D  is  the  path  length.  This  expression  gives  the  predicted  field 
strength  in  these  cases. 

As  already  mentioned,  the  assumptions  of  geometric  ray  propagation  are 
not  valid  for  long  path  lengths.  At  these  distances  the  easiest  procedure  is 
to  use  an  empirical  model  which,  hopefully,  has  been  constructed  from  a  data 
base  containing  examples  of  all  the  loss  mechanisms  accounted  for  in  the 
short-path  model  plus  those  losses  unique  to  long-path  propagation.  The  only 
way  we  can  ascertain  whether  the  model  contains  all  the  loss  mechanisms  is  by 
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comparing  predictions  against  actual  measured  data.  The  results  of  such  a 
comparison  appear  later  in  this  report,  but,  first,  we  describe  the  empirical 
model  tested. 

Path  Lengths  Greater  Than  9000  km  (the  FTZ  Model).  The  FTZ  method  of 
field  strength  calculation  was  developed  by  Beckmann  (1967)  and  has  been 
continually  refined  by  the  work  of  Damboldt  (1975)  at  FTZ.  At  these  path 
lengths,  the  effects  of  ionospheric  tilts  and  scattering  from  irregularities 
cause  contributions  to  the  measured  field  strength  from  paths  outside  the  main 
great  circle  path.  These  effects  also  complicate  assigning  an  antenna  gain  to 
a  particular  path  since  energy  may  be  received  from  many  azimuths.  These 
considerations  have  led  to  the  development  of  the  following  empirical  model, 
which  is  used  by  the  CCIR  for  paths  greater  than  9000  km  and  is  used,  in  an 
altered  form  to  be  described  later,  for  all  path  lengths  by  PROPHET  programs. 

Another  propagation  mode  that  is  always  present  to  some  degree  is  at 
frequencies  exceeding  the  MUF  for  the  path.  The  ray  path  geometry  methods 
described  earlier  would  predict  that  a  ray  at  such  a  frequency  would  penetrate 
the  ionosphere  and,  thus,  not  be  received  by  a  ground-based  receiver.  It  has 
been  found  experimentally  however  that  usable  field  strength  is  available  at 
frequencies  that  exceed  the  MUF.  Some  factors  that  contribute  to  this 
phenomenon  are  forward  scatter  via  F  reflection,  lateral  ground  scatter  from 
great  circle  paths  with  a  higher  F2-layer  MUF,  sporadic-E  scatter,  supermodes, 
and  extremely  long  Pederson  (high  angle)  rays  (Beckman,  1967). 

Beckmann  determined  that  the  measured  field  strength  was  a  function  of 
both  the  lowest  usable  frequency  (fj^)  for  the  path,  which  he  defines  as  that 
frequency  for  which  the  field  strength  is  some  low  limiting  value,  and  the 
MUF.  Experiments  showed  that  the  field  strength  rises  with  frequency  from  f^ 
to  a  peak  at  some  intermediate  frequency  somewhat  below  the  MUF  and  then 
decreases  to  the  same  low  limiting  value  at  a  frequency  somewhat  greater  than 
the  MUF,  which  he  calls  the  operational  MUF.  The  operational  MUF  (f^)  is 
obtained  as  the  product  of  a  numerical  factor,  determined  from  geographical 
and  diurnal  characteristics  of  the  path,  times  the  MUF.  The  f^^  -  f^ 
combination  then  defines  the  transmission  frequency  range  (TFR)  for  the  path, 
or  that  range  of  frequencies  for  which  usable  field  strength  may  exist, 
depending  on  system  parameters  such  as  noise  level  and  receiver  sensitivity. 

a.  FTZ  Long- Path  Model.  Given  the  above  behavior  of  the  field  strength 
with  frequency,  Beckmann  (1967)  used  the  following  equation 
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to  represent  the  rise  of  field  strength,  with  frequency,  from  f  .  Here  F  is 

o 

the  free-space  field  strength  given  by 


F 

o 


20  log 


3  X  10  ./Perp 


(17) 


where  Perp  is  the  effective  radiated  power  in  kilowatts  and  d  is  the  slant 
path  range  in  kilometers  for  mirror  reflection  from  a  height  of  300  km. 

Equation  16  is  exactly  valid  for  nondeviative  absorption  and  can  be 
derived  in  the  following  way.  As  we  saw  earlier,  the  D- region  (nondeviative) 
absorption  is  inversely  proportional  to  the  square  of  the  operating  frequency. 
At  low  frequencies,  the  D-region  absorption  is  greatest  and  is  the  dominant 
factor  in  determining  operational  frequencies.  So  we  can  approximate  the 
field  strength  at  the  low  frequencies  by 


where  is  the  ionospheric  absorption  discussed  earlier. 
11,  has  the  approximate  form  for  oblique  propagation 


(18) 

Now  from  equation 
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(f  *  f„)'  ' 


(19) 


where  we  have  replaced  the  longitudinal  component  of  the  gyrofrequency  with 
the  full  gyrofrequency,  which  introduces  only  a  slight  error.  We  then  define 

3 

as  the  frequency  that  produces  a  field  strength  of  0  dB  (;uV/m)  for  10  kW 
effective  radiated  power.  Equation  18  then  becomes 


0  =  F  (10  kW)  - 
o 


A  sec  i 


(20) 


or 


F  (10^  kW)(f^  +  f  )^  =  A  sec  i  . 
o  L  H  o 

Putting  this  expression  into  equation  19,  equation  18  becomes 

2, 


(21) 
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Although  exact  for  nondeviative  absorption,  use  of  this  expression  in 
general  is  adequate  since  at  low  frequencies  nondeviative  absorption  dominates 
other  loss  mechanisms. 

To  include  the  effects  of  the  operational  MUF  in  the  field  strength, 
Beckmann  (1967)  introduced  the  following  changes  in  equation  16  to  produce  the 
decline  in  field  strength  with  frequency  from  the  peak 


^ti  =  ^o 


[<£  *  +  fH>'l 

1  '  c 

(22) 


where 


Note  this  expression  for  field  strength  satisfies  the  requirement  that  for  f  = 

3 

and  f  =  f^  the  field  strength  for  10  kW  effective  radiated  power  is  0  dB. 

The  field  strength  peaks  at  a  frequency  fp  given  by  the  geometric  mean  of  f^^ 

and  f  , 
o 

fp  =  y(fj^  +  f^)(f^  +  fj^)  . 

To  derive  the  expression  for  the  field  strength  from  f  to  the  peak,  we 
assumed  an  effective  radiated  power  of  1000  kW.  Equation  22  can  be  adjusted 
for  arbitrary  power  levels  and  antenna  gains  by  including  the  actual  power 
level  used  and  by  addition  of  the  transmitter  antenna  gain.  For  example,  for 
100  kW  radiated  power  we  have 

r(f  +  f  )^  (f  +  f  )^'j 

=  (139.6  -  20  log  d)  1  -  C  - ^  ^  -  14.8  + 

[  l(fo  +  V  ^  J  J 

where  is  the  maximum  transmitting  antenna  gain  in  dBi  in  the  range  of  8°  to 
10°  elevation  angle.  Since  Beckmann  originally  used  a  reference  field  due  to 
an  isotropic  antenna  over  perfect  ground,  which  has  a  gain  of  4.8  dB  with 
respect  to  an  isotropic  antenna  in  free  space,  4.8  is  subtracted  to  go  from 
effective  radiated  power  to  radiated  power. 

Likewise,  the  field  strength  for  a  1-kW  transmitter  would  be 


f<£  +  <f,  +  f„)2l 

1  -  c 

-  2  2 

=  (139.6  -  20  log  d) 


In  general,  for  arbitrary  radiated  power,  we  have 
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=  (139.6  -  20  log  d) 
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1000 


(23) 


where  P,^  is  the  transmitted  power  in  watts. 

The  operational  MUF  used  by  Beckman  has  been  changed  to  include  more 
geographical  and  diurnal  effects  by  Damboldt  (1975).  The  constants  that 
appear  in  the  expression  below  have  been  changed,  and  continue  to  change,  as 
comparison  with  data  warrants.  For  completeness,  we  also  describe  the  FTZ 
extension  to  paths  less  than  7000  km,  although  this  is  not  used  in  the  HFBC84 
program. 

The  current  expression  for  the  operational  MUF  is  (Bradley  and  Liu,  1982) 

f  =  K  •  MUF 
o 

where,  for  E-layer  modes,  K  =  1.  For  F2-layer  modes 
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with  K  =  1  for  total  path  ground  range,  D,  greater  than  3600  km  and 


(24) 


K  =  1.25  +  .25  cos 
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for  D  <  3600  km. 
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In  equation  24,  f  ,  f  ,  and  f  ,  are,  respectively,  the  MUF  for  the 
g  g,max  g.rain 

hour  of  interest  and  the  maximum  and  minimum  MUF  for  the  24 -hour  day.  W,  X, 
and  Y  are  empirical  constants,  which  are  functions  of  raypath  azimuth  at  the 
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midpath  point  derived  by  linear  interpolation  with  angle  between  W  -  .1,  X  - 
1.2,  Y  =  6  for  east-west  paths  and  W  -  .2,  X  -  .2,  Y  =  .4  for  north-south 
paths . 

To  determine  the  operational  MUF  for  a  particular  path,  the  general 
technique  is  as  follows:  first,  the  ionospheric  reflection  positions  of  the 

least  order,  equal  hop  length  F  mode  having  a  hop  length  of  less  than  3600  km 

are  determined,  along  with  those  of  the  least  order,  equal  hop  length  E-mode 
with  a  hop  length  of  less  than  1800  km.  If  the  ground  range  will  support  the 
single -hop  F  mode,  the  ionospheric  reflection  is  taken  at  the  midpath  and  K  is 
determined  at  that  point.  Otherwise,  four  control  points  are  defined,  each 
being  a  position  of  ionospheric  reflection  for  the  E  and  F  modes  closest  to 
the  transmitter  and  receiver.  At  each  control  point  the  K  factor  is 
determined  by  the  above  method,  and  the  operational  MUF  for  the  hop  is 
determined  at  each  control  point.  Then  the  larger  of  the  two  operational  MUFs 
is  chosen  as  the  value  for  the  path.  The  minimum  of  these  two  values  is  used 

as  the  value  of  f^  for  the  path  at  that  hour  (Bradley  and  Liu,  1982) . 

The  other  input  to  the  field  strength  expression  is  the  frequency  f^. 
Referring  to  equation  20  we  have 


A  sec  6 
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and  we  see  that  the  fj^  depends  on  the  ionospheric  absorption  loss  term.  Since 
Beckmann  originally  developed  the  technique,  the  actual  equation  used  has 
changed  several  times,  again  as  new  developments  have  demanded  (Damboldt, 
1975).  The  form  of  the  equation  is,  however,  still  that  given  in  equation  19. 
Currently,  the  expression  used  for  the  daytime  f^  in  the  FTZ  models,  for  path 
lengths  greater  than  2000  km  (Bradley  and  Liu,  1982),  is 


(5.3  .  I) 


'  1  + 

.009R 

12 

20  log[ 

3.5  • 

d 

Icos 

1/2 


Z  ycosx. 

i=l 


1/2 


W 


(25) 


where  for  F  and  E  modes  alike: 

I  “  seasonal  factor  with  values  in  the  range  of  1.0  to  1.1,  depending  on 
the  transmitter/receiver  latitude  combination  and  month 
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R^2  “  12 -month  running  mean  sunspot  number 

4)  =  angle  of  incidence  at  90  km,  assuming  mirror  reflection  from  300  km 

and  hop  span  D  for  the  great  circle  path  ray 
d  =  slant  path  length  for  mirror  reflection  from  300  km 
N  =  niimber  of  legs  of  least  order  mode  (twice  the  number  of  hops) 

=  solar  zenith  angle  at  a  height  of  90  km  where  the  ith  leg  of  the 

raypath  achieves  a  height  of  90  km 
W  =  winter  anomaly  term  determined  at  path  midpoint 

W  is  unity  for  geographic  latitudes  of  0°  to  30°  and  at  90° ,  and  peaks  at 
60°  latitude  with  values  ranging  from  1.0  to  1.30  depending  on  season.  If 
some  of  the  hop  legs  are  in  darkness,  x^  is  taken  as  90°  in  the  summation  for 
those  legs  (Bradley  and  Liu,  1982). 

At  night  the  solar -driven  ionization  process  ceases  and  the  absorption  is 
greatly  reduced.  During  this  period,  the  frequency  fj^  is  given  by 
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73000 


where  D  is  ground  range  in  km.  For  the  evening  transition  hours,  f^^  is 
assumed  to  fall  exponentially  with  time  from  a  value  f^  =  2  f^  at  a  reference 
time  T  to  twice  the  nighttime  value,  as 


f^  =  2  f^  exp[-.23(t  -  T)] 

giving  a  decay  to  f^^^  in  roughly  3  hours  (Bradley  and  Liu,  1982). 

This  completes  the  specification  of  the  FTZ  long-path  model.  As 
mentioned  above,  the  CCIR  recommends  use  of  this  method  exclusively  for  paths 
greater  than  7000  km.  At  these  long  distances,  we  are  also  liable  to  see 
focusing  of  energy  due  to  layer  tilts  and  irregularities.  This  results  in  an 
increase  of  measured  field  strength.  This  effect  is  included  in  HFBC84  with 
the  addition  of  a  term,  G  ,  to  equation  23.  For  propagation  to  distances  D 
(great  circle  distance)  greater  than  7rR/2,  where  R  is  the  earth's  radius, 
focusing  is  taken  into  account  through  the  formula 


%  =  -  20  log^l 
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[-^2  >  D  >  ~]  »rR,  n  =■  1  and  2 
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As  tends  to  infinity  for  D  =  njrR,  its  value  is  limited  to  30  dB  {Report  to 

the  Second  Session,  1984). 

FTZ  originally  used  the  model  we  have  outlined  for  all  path  lengths,  and 
all  PROPHET  programs  now  use  a  modified  version  of  the  model.  Recently, 
however,  FTZ  has  developed  an  improved  model  for  paths  less  than  2000  km  and, 
together  with  the  above  long-path  model,  covers  the  entire  range  of  path 
lengths.  Although  this  model  is  not  implemented  in  HFBC84  or  PROPHET,  we 
present  it  here  for  reference  and  as  a  possible  addition  to  future  PROPHET 
field  strength  prediction  models. 

b.  FTZ  Short-Path  Model.  The  form  of  the  model  for  short  paths 
(<2000  km)  is  similar  to  the  long-path  model  with  the  following  modifications. 
An  important  input  is  the  frequency  f„.„  at  which  the  field  strength  is  a 

MAJv 

maximum.  This  frequency  is  given  by  (Bradley  and  Liu,  1982) 
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where  f  is  the  MUF  for  the  1-hop  E  or  F2  mode  to  ground  range  D,  whichever  is 
g 

greater.  Note  that  D/3000  is  the  square  of  the  nighttime  f^. 

The  operational  MUF  is  given  by 
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and  again  we  have  different  equations  for  f  depending  on  time  of  day. 

For  daytime  periods ,  the  f is  given  by 

fjjj  =  1.1  •  yi  -  0.003  •  R12  ♦  fj^ 

where  f^  is  the  daytime  fj^  used  for  the  long-path  model  assuming  a  1-hop  F2 
mode  reflected  from  a  height  of  300  km.  For  nighttime  and  evening  transition 
hours,  the  methods  outlined  for  the  long-path  model  are  used  with  the  above 
expression  for  the  daytime  frequency  f^^,  replacing  the  earlier  expression  in 
the  calculation  of  the  evening  transition  f  . 

The  field  strength  for  f  <  f„.„  is  then  given  by 
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and  for  f  >  f,,.„  we  have 
MAX 


The  parameters  and  are  the  same  as  in  the  long-path  model. 

Over- the -MUF  Field  Strength  Model.  As  mentioned  earlier,  the  observed 
field  strength  does  not  abruptly  vanish  as  the  operating  frequency  begins  to 
exceed  the  basic  MUF  for  a  path.  In  fact,  measurable  field  strengths  are 
frequently  observed  at  frequencies  two  or  more  times  the  MUF. 

If  we  ignore  MUF  prediction  errors  and  the  monthly  median  nature  of  the 
parameters  involved,  we  can  identify  several  factors  that  can  contribute  to 
this  over -the -MUF  propagation.  Some  of  these  mechanisms  are  blanketing 
sporadic  E  layers  that  support  skywave  propagation  at  higher  frequencies, 
ground  scatter  from  off  great  circle  paths  that  may  have  larger  MUFs,  and 
forward  scatter  via  F-layer  reflection  (Beckmann,  1967). 

As  the  operating  frequency  rises,  some  field  strength  reduction  is 
noticeable  even  at  frequencies  below  the  path  MUF.  This  is  because 
frequencies  near  the  MUF  penetrate  deeper  into  the  reflecting  layer  and  are 
subject  to  more  scattering  out  of  the  great  circle  path.  The  peak  field 
strength  is  then  seen  at  a  frequency  somewhat  below  the  MUF. 

There  is  experimental  evidence  that  the  rate  of  decrease  of  field 
strength  at  frequencies  greater  than  the  MUF  may  also  be  dependent  on  path 
length  (Dieminger  and  Rose,  1961).  This  is  because  for  a  longer  path  length  a 
wave  is  more  likely  to  encounter  ionization  density  irregularities  and 
ionospheric  tilts.  So  there  is  an  increased  capability  for  energy  at  higher 
frequencies  to  be  scattered  to  the  receiver.  This  results  in  measurable  field 
strength  at  increased  frequencies  and  a  slower  roll-off  with  frequency  from 
the  peak  field  strength. 
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Overall  then,  we  can  have  usable  field  strength  at  frequencies  greater 
than  the  path  MUF,  although  the  level  of  the  field  strength  is  highly 
variable.  This  variability  depends  on  the  state  of  the  ionosphere,  the  length 
of  the  path,  and,  possibly,  on  terrain  in  the  case  of  ground  scatter. 

Some  of  these  considerations  are  inherent  in  the  empirical  FTZ  method  of 
field  strength  prediction.  By  employing  the  operational  MUF,  recognition  is 
given  to  the  fact  that  there  is  propagation  at  these  higher  frequencies.  Some 
controversy  remains,  however,  as  to  the  correct  form  this  roll-off  in  field 
strength  takes  in  this  model.  We  shall  have  more  to  say  on  this  subject  in  the 
discussion  section. 

In  our  earlier  discussion  of  the  HFBC84  model  for  paths  less  than  7000 
km,  we  did  not  discuss  the  over- the -MUF  loss  turn  in  equation  3.  We  now  take 
up  this  topic  and  present  the  model  used  in  that  program.  This  model,  as  are 
most  over-the-MUF  loss  calculations,  is  based  on  tbe  Philips-Abel  theory 
(Wheeler,  1966). 

In  this  theory,  the  reflecting  layers  are  considered  to  be  composed  of 
individual  ionization  patches.  Each  patch  consists  of  subpatches  with  varying 
ionization  levels,  for  each  of  which  there  exists  a  classical  MUF.  The  median 
of  the  MUFs  for  the  subpatches  is  the  MUF  for  a  patch.  The  number  of  patches 
supporting  wave  reflection  falls  with  increasing  frequency  and  so  there  is  no 
abrupt  signal  cutoff  (Wheeler,  1966). 

Wheeler  (1966)  has  suggested  that  for  a  wave  of  frequency  f  greater  than 
the  MUF,  the  median  signal  power  can  be  taken  as  proportional  to  the  number  of 
patches  of  F- region  ionization  that  support  transmission  at  f.  That  is, 
proportional  to  the  number  of  patches  such  that  f<MUF.  With  the  frequency 
distribution  of  the  MUFs  of  the  sub-patches  taken  to  be  Gaussian  we  have 
(Report  by  the  Chairman,  1975) 
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and  X  = 
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Here  o  is  the  standard  deviation  of  the  MUF  distribution 


o  =  .78  (F^  -  1)  •  MUF 
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where  is  the  90-percentile  (highest  possible  frequency  (HPF))  factor  for 
the  MUF  distribution  (Barghausen  et  al.,  1969),  and  the  MUF  is  the  median 
value  predicted  by  the  methods  presented  earlier. 

A  much  simplified  method  based  on  the  model  described  above  is  used  in 
HFBC84.  This  method  is  based  on  a  fit  to  a  limited  set  of  over- the -MUF  data 
collected  by  Wheeler  (Bradley  et  al . ,  1982).  The  loss  term  is  given  by 


and  is  not  allowed  to  become  larger  than  85  dB.  There  is  some  question  as  to 
the  accuracy  of  this  simple  method,  and  we  will  have  more  to  say  on  this 
subject  in  the  discussion  section. 

This  completes  the  specification  of  the  HFBC84  field  strength  model.  We 
have  given  a  detailed  description  to  point  out  the  various  mechanisms  for 
energy  reduction  that  contribute  to  the  field  strength. 

The  parent  program  of  HFBC84  contains  a  more  elaborate  ionospheric  loss 

calculation  and  also  contains  several  more  loss  terms  in  the  field  strength 

equation  (polarization  loss,  E  obscuration  loss,  etc.)  (Second  CCIR  Computer- 

s 

based,  1978).  Since  they  are  not  included  in  HFBC84,  we  will  not  discuss  them 
here.  The  reader  wishing  information  about  these  loss  mechanisms  is  directed 
to  the  aforementioned  reference. 

Next  we  describe  the  field  strength  model  used  in  the  PROPHET  family  of 
propagation  prediction  programs.  This  model  is  based  on  the  FTZ  method,  which 
we  have  already  described.  In  this  case,  the  FTZ  model,  or  the  modified 
version  of  it  used  in  PROPHET,  is  used  at  all  path  lengths.  Further,  all 
references  to  the  MUF  will  be  understood  to  refer  to  that  value  predicted  by 
the  MINIMUF  algorithm.  Since  there  is  currently  no  E- layer  model  in  PROPHET, 
these  values  will  refer  to  F2- layer  MUPs  only. 

PROPHET  Field  Strength  Model 

The  model  employed  here  is  the  same  as  that  presented  for  the  long-path 
method  used  in  HFBC84,  with  several  modifications. 

The  model  employed  in  PROPHET  uses  the  ionospheric  absorption  index  used 
in  the  QLOF  routine  in  PROPHET.  This  is  an  empirical  LUF  prediction  program 
and  uses  an  ionospheric  absorption  term  based  on  the  work  of  Schultz  and 
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Gallet  (1970)  to  predict  that  parameter.  The  expression  used  in  QLOF,  which 
replaces  the  index  Ai  in  equation  11,  is 


A. 
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CH  (921.0.V)  -|-2-Q  •  M 

CH  (921.0,  -  <!>)] 


where 

A  =  286  •  W  .  (1.0  +  .5U|)  cos^(^S  -  i)  .  (29) 

s 

Here  CH  is  the  Chapman  function,  is  the  subsolar  geographic  latitude, 

(f>  is  geographic  latitude,  and  x  is  the  solar  zenith  angle  measured  at  W  is 
a  winter  anomaly  term  defined  by 


W  -  1.0  +  0.0275[30.0  -  | (60  -  |^1)|]  , 
M  is  defined  by. 
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For  paths  less  than  4000  km,  all  calculations  are  done  at  a  single 
control  point,  the  path  midpoint.  For  paths  greater  than  4000  km,  control 
points  are  established  at  the  midpoint  and  at  points  1000  km  from  each  end 
point.  A^  is  then  determined  at  each  control  point,  and  the  average  of  these 
values  determines  the  absorption  index  for  the  path.  A  more  detailed 
description  of  QLOF  is  available  in  Argo  and  Hill  (1977)  and  Sailors  and 
Moision  (1987). 

The  LUF-MUF  combination  defines  the  rise  in  field  strength  with  frequency 
up  to  the  MUF,  as  outlined  earlier.  Now  we  must  determine  the  operational 
MUF,  which  appears  in  the  FTZ  model  and  relates  to  the  field  strength  we 
observe  at  frequencies  larger  than  the  MUF. 
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As  we  have  discussed  above,  the  field  strength  measured  at  frequencies 
greater  than  the  MUF  is  due  to  various  scattering  mechanisms ,  both  in  the 
ionosphere  and  on  the  earth's  surface.  As  such,  the  Wheeler  method  of 
predicting  over-the-MUF  loss,  which  we  described  earlier,  bears  no  physical 
connection  to  the  true  causes  of  this  type  of  propagation.  It  does,  however, 
provide  a  convenient  means  of  treating  this  phenomenon  in  terms  of  the 
relevant  parameter,  i.e.,  frequency  to  MUF  ratio. 

The  data  on  which  the  empirical  PROPHET  field  strength  model  and  the  FTZ 
model  are  based  were  collected  over  a  number  of  years  at  the  FTZ  facility. 
Since  there  is  likely  to  be  geographic  dependence  in  over-the-MUF  propagation 
phenomenon,  we  believe  the  operational  MUF  determination  used  in  the  FTZ  model 
may  contain  aspects  that  are  specific  to  that  site. 

For  these  reasons,  a  simplified  method  of  determining  the  operational  MUF 
has  been  implemented  in  the  PROPHET  field  strength  calculation.  As  in  the 
Wheeler  method,  we  assume  that  the  MUF  value  at  a  given  time,  season,  geo¬ 
graphic  position,  and  sunspot  number  is  a  statistical  parameter  and  has  a 
Gaussian  distribution  around  a  mean.  This  mean  value  is  that  predicted  by 
MINIMUF.  Given  this  assumption,  the  operational  MUF  in  PROPHET  is  determined 
by  the  product  of  the  MUF  and  the  99.1  percentile  value  of  the  MUF 
distribution 

f  -  1.85  •  (HPF)  ♦  MUF, 
o 

where  HPF  is  the  90-percentile  value  factor  for  an  assumed  Gaussian  MUF 
distribution  (Barghausen  et  al.,  1969).  This  technique  for  determining  f^ 
provides  a  value  that  is  somewhat  between  the  two  methods  outlined  earlier  in 
its  philosophy.  Its  adequacy  for  determining  f^  and,  hence,  for  predicting 
field  strength  at  frequencies  greater  than  the  MUF  will  be  discussed  in  the 
data  comparison  results  presented  next. 

With  these  differences,  the  field  strength  calculation  in  PROPHET  is  the 
same  as  that  described  for  the  long-path  model  in  HFBC84. 
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FIELD  STRENGTH  DA'^A  BASE  AND  SCREENING  PROCEDURE 


For  the  comparison  of  predicted  values  to  measured  data,  we  employed  a 
statistical  data  screening  program,  DASCR3 ,  which  has  been  described  in 
several  previous  NOSC  technical  reports  (Sailors  et  al.,  1981  and  1986).  This 
program  allows  the  data  base  to  be  separated  at  run  time  into  blocks  that 
satisfy  certain  screening  conditions.  In  this  way  we  can  delineate  subsets  of 
a  large  data  base  for  independent  testing.  Using  this  feature  of  the  program, 
we  partitioned  the  data  base  into  blocks  corresponding  to  path  length  ranges, 
sunspot  ranges,  time  ranges,  etc.  These  were  then  screened  separately  to  test 
accuracy  under  these  conditions. 

The  data  base  used  in  this  analysis  is  a  modified  version  of  CCIR  Data 

uV 

Base  C.  It  consists  of  12,277  median  skywave  field  strength  (dB  above  1 
data  values  normalized  to  1  kW  erp.  This  differs  from  the  full  Data  Base  C 
in  that  all  "long  paths"  have  been  removed.  These  are  paths  for  which  the 
radiowave  travels  the  long  way  around  the  world.  This  was  done  because 
PROPHET  does  not  currently  have  the  option  to  calculate  field  strength  on  that 
type  of  path. 

Table  2  gives  a  complete  summary  of  the  measured  data  set  used  in  this 
analysis.  Shown  here  are  circuit  identification  number,  transmitter  and 
receiver  site  names  and  geographic  coordinates,  midpath  latitude,  azimuth  from 
transmitter  to  receiver,  path  length,  number  of  path  hours,  and  frequencies 
for  each  circuit.  We  include  circuit  number  and  midpath  latitude  for  later 
referral,  as  these  are  variables  used  as  a  means  to  partition  the  data  base 
for  screening. 

The  data,  collected  from  1964  to  1984,  represent  almost  two  complete 
solar  cycles  and  contain  data  samples  representing  12 -month  running  mean 
sunspot  numbers  in  the  range  13  to  165. 

For  each  circuit  in  table  2  there  are,  in  general,  several  frequencies 
for  which  data  were  obtained,  and  these  are  also  indicated.  In  some 
instances,  the  frequency  was  greater  than  the  MUF  determined  for  that  circuit¬ 
time  combination.  Thus  we  could,  to  some  extent,  separately  investigate  the 
over-the-MUF  models  used  by  the  respective  programs.  Since  there  are 
differences  in  the  HUF  prediction  routines,  a  particular  frequency  may  be  over 
the  MUF  for  one  model  and  under  the  MUF  for  the  other.  Thus,  population  sizes 
will  differ  between  models  in  these  comparison  results. 
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Table  2.  Summary  of  transmission  circuits  contained  in  the  field  strength  data  base. 
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Table  2.  Summary  of  transmission  circuits  contained  in  the  field  strength  data  base 
(continued) . 
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The  data  base  used  to  assess  these  models,  Data  Base  C  from  the  CCIR,  has 
no  frequencies  below  4.8  MHz.  There  is  currently  a  new  data  base  available 
from  the  CCIR,  called  Data  Base  D,  which  has  frequencies  down  to  about  2  Hz. 

We  recommend  that  this  data  base  be  obtained  for  further  efforts . 

The  procedure  followed  was  to  generate  separate  predicted  data  bases 
using  each  of  the  models  described  in  this  report.  These  were  then  compared, 
using  DASCR3,  to  the  measured  data  in  blocks  corresponding  to  specified 
screening  conditions  as  described  earlier. 
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RESULTS  OF  COMPARISON  TO  DATA 


We  present  here  the  results  of  the  comparison  of  the  PROPHET  and  HFBC84 
field  strength  prediction  programs  to  measured  data.  Throughout  this  section, 
table  2  should  be  consulted  as  the  relative  lack  of  data  in  some  cases  makes 
the  validity  of  statistical  analysis  questionable.  A  limited  analysis  of  the 
results  and  recommendations  for  further  investigations  will  be  presented  in 
the  next  section. 

The  following  quantities  are  referred  to  in  the  accompanying  figures 
(Sailors  et  al . ,  1981); 

•  residual  =  observed  datum  -  predicted  value 

•  relative  residual  =  residual/observed  datum 

•  absolute  relative  residual  “’absolute  residual/observed  datum. 

Here  the  observed  datum  and  predicted  value  correspond  to  a  complete  set 
of  identical  conditions,  i.e.,  time,  month,  year,  sunspot  number,  path,  etc. 
From  these  parameters,  the  following  six  statistical  quantities  derived  from 
the  residual  distributions  are  plotted: 

•  average  residual  (bias) 

•  root  mean  square  residual  (rms  res) 

•  average  relative  residual  (av  rel  res) 

•  root  mean  square  relative  residual  (rms  rel  res) 

•  average  absolute  relative  residual  (av  abs  rel  res) 

•  correlation  coefficient  between  observed  and  predicted  values. 

Each  of  these  statistical  quantities  is  calculated  by  the  DASCR3  program. 
Examples  are  shown  in  figures  1  to  6  where  each  of  these  quantities  is  plotted 
as  a  function  of  the  operating  frequency  to  MUF  ratio. 

OVERALL  RESULTS  AS  A  FUNCTION  OF  F/MUF  RATIO 

As  we  pointed  out  earlier,  there  is  some  question  as  to  the  adequacy  of 
current  models  for  field  strength  prediction  when  the  frequency  exceeds  the 
MUF  for  the  mode.  The  series  of  plots  in  figures  1  to  6  show  clearly  that  for 
f/MUF  ratios  of  -1.5  or  larger  neither  model  provides  satisfactory  predictions. 

The  average  residual,  figure  1,  and  the  average  relative  residual,  figure 
2,  locate  the  center  of  the  distributions  for  these  quantities  and  are 
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Figure  4.  Root  mean  square  residual  as  a  function  of  frequency-to-MUF  ratio. 
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referred  tc  as  the  bias  in  the  estimate  (Sailors  et  al.,  1981).  Both  of  these 
figures  show  large  errors  at  ratios  greater  than  approximately  1.5.  Also,  the 
PROPHET  model  shows  a  relative  peak  in  the  error  at  low  frequencies. 

The  average  absolute  relative  residual  is  a  measure  of  the  average 
magnitude  of  the  error  (bias) .  Figure  3  shows  the  same  characteristic 
behavior  as  the  earlier  figures.  Both  models  show  quite  good  results  for 
frequencies  of  approximately  0.5  to  1.5  of  the  MUF.  However,  both  models  show 
poor  performance  at  higher  frequencies.  Again,  the  PROPHET  model  shows  some 
problems  at  low  frequencies. 

The  root  mean  square  residual  and  root  mean  square  relative  residual  are 
measures  of  the  dispersion  in  the  error  (Sailors  et  al.,  1981).  They 
represent  the  standard  deviation  of  the  error  about  the  origin.  So  for  a 
model  with  zero  bias  in  the  error,  the  standard  deviation  and  rms  values  are 
identical.  Figures  4  and  5  show  these  quantities  as  a  function  of  f/MUF 
ratio.  They  both  display  large  dispersion  in  the  predictions. 

The  correlation  coefficient  is  a  measure  of  the  degree  of  association  of 
the  closeness  of  fit  between  the  observed  and  predicted  values.  It  indicates 
the  strength  of  the  tendency  for  high  (low)  values  of  one  variable  to  be 
associated  with  high  (low)  values  of  the  other  (Sailors  et  al.,  1981).  Figure 
6  shows  the  correlation  coefficient  in  the  current  case.  Both  models  show 
similar  behavior  in  that  they  tend  to  decrease  almost  uniformly  as  the 
frequency  to  MUF  ratio  increases.  HFBC84  shows  an  offset  that  increases  its 
correlation  over  most  of  the  frequency  range,  but  it  shows  the  same  overall 
tendency  as  PROPHET . 

These  results  show  that,  in  general,  neither  method  of  predicting  over- 
the-MUF  field  strength  is  adequate.  For  this  reason  the  following  plots  and 
tables  have  been  screened  according  to  whether  the  operating  frequency  is 
greater  or  less  than  the  MUF.  Results  are  presented  for  each  case.  In  this 
way  other  aspects  of  the  models'  performance  can  be  investigated,  independent 
of  the  effects  of  the  poor  over-the-MUF  models.  In  the  appendix  we  give  the 
overall  results  in  tabular  form. 
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INDIVIDUAL  CIRCUIT  RESULTS 

In  figures  7  to  18  we  show  the  performance  of  the  models  for  the 
individual  circuits  listed  in  table  2.  In  figures  7  to  12  the  operating 
frequencies  are  less  than  (or  equal  to)  the  MUF  for  the  path,  and  in  figures 
13  to  18  the  frequencies  are  greater  than  the  MUF. 

Figure  7  shows  considerable  variation  (+20  dB)  in  the  average  residual 
(bias)  around  a  value  near  zero  for  most  circuits.  However,  for  circuit 
identification  numbers  less  than  about  18,  PROPHET  tends  to  underpredict  the 
field  strength  relative  to  HFBC84  predictions.  For  circuit  numbers  55  and 
grater,  HFBC84  tends  to  underpredict  relative  to  PROPHET.  Referring  to  table 
2,  circuit  numbers  less  than  18  refer  to  paths  less  than  1200  km,  while  those 
greater  than  55  refer  to  paths  approximately  7000  km  and  longer.  Since  PROPHET 
uses  a  variation  of  the  FTZ  method  at  all  path  lengths,  both  of  these  path 
length  regions  point  to  a  problem  with  that  model. 

The  average  residual  for  frequencies  greater  than  the  MUF  is  shown  in 
figure  13.  The  tendency  here  is  for  both  models  to  underpredict  at  short 
(<2000  km)  path  lengths .  Both  models  overpredict  for  longer  (>4000  km)  paths 
with  PROPHET  showing  a  decreasing  linear  trend  with  path  length. 

In  figures  8  to  12,  the  plots  for  the  remaining  statistical  parameters 
are  presented  for  frequencies  less  than  the  MUF.  Figure  8  shows  that  the  rras 
residual  averages  approximately  10  dB  for  short  to  midlength  paths  for  both 
models.  Again  we  see  a  small  long  path  deviation  in  HFBC84. 

Figures  9  to  11  all  show  similar  behavior.  The  dominant  feature  is  the 
erratic  performance  for  circuit  22.  Although  the  data  sample  is  extremely 
small  for  this  circuit,  the  failure  of  both  models  tends  to  make  the  quality 
of  the  data  suspect. 

Figure  12  shows  the  correlation  coefficient  for  each  model.  Both  models 
show  variation  over  the  entire  range  from  1.0  to  -1.0.  No  clear  trend  is 
discernible  from  this  plot. 

For  completeness,  figures  13  to  18  show  the  results  for  frequencies 
greater  than  the  MUF.  As  we  have  indicated,  the  poor  performance  of  the  over- 
the-MUF  models  in  general  make  isolating  other  effects  difficult. 
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Figure  11.  Average  absolute  relative  residual  for  each  circuit  with 
frequency  less  than  or  equal  to  the  MUF. 
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Figure  12.  Correlation  coefficient  for  each  circuit  with  frequency  less 
than  or  equal  to  the  MUF. 
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Figure  13.  Average  residual  for  each  circuit  with  frequency  greater  than 
the  MUF. 
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Figure  14.  Root  mean  square  residual  for  each  circuit  with  frequency 
greater  than  the  MUF. 
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Figure  16.  Root  mean  square  relative  residual  for  each  circuit  with 
frequency  greater  than  the  MUF. 
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Figure  17.  Average  absolute  relative  residual  for  each  circuit  with 
frequency  greater  than  the  MUF. 
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Figure  18.  Correlation  coefficient  for  each  circuit  with  frequency  greater 
than  the  MUF. 
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CIRCUIT  LENGTH  RESULTS 


The  results  presented  in  the  previous  section  indicate  the  relatively 
poor  performance  of  the  HFBC84  model  for  paths  longer  that  -7000  km.  The 
results  presented  here  are  consistent  in  that  HFBC84  again  shows,  in  figure 
19,  an  increased  bias  for  these  path  lengths  at  frequencies  below  the  MUF. 
Except  for  the  very  longest  paths  (>11,000  km),  PROPHET  is  relatively 
consistent  with  an  average  bias  of  +5  dB  across  all  range. 

Results  for  the  rms  residual  are  shown  in  figure  20.  Both  models  vary 
between  approximately  10  to  20  dB  over  roost  of  the  ranges.  PROPHET  shows  even 
larger  variation  at  the  extreme  path  lengths. 

Figure  21  shows  the  average  relative  residual  for  frequencies  less  than 
the  MUF.  Again  we  see  HFBC84  performs  relatively  poorly  at  longer  path  ranges. 
In  this  case,  PROPHET  also  shows  large  variation  at  longer  path  lengths. 

Results  for  the  other  statistical  parameters  for  frequencies  less  than 
the  MUF  are  shown  in  figures  22  to  24.  In  particular,  the  correlation 
coefficient,  figure  24,  shows  a  large  amount  of  variation  for  both  models.  No 
clear  systematic  trends  are  discernible  from  these  results. 

Figures  25  to  30  show  plots  of  the  statistical  parameters  for  frequencies 
greater  than  the  MUF.  The  broken  lines  in  these  figures  indicate  a  lack  of 
data  for  that  particular  model.  Neither  model  performs  well  in  this  series  of 
plots.  One  feature  that  does  stand  out  is  the  similarity  in  the  results  in 
each  case,  which  makes  specifying  the  "better”  model  for  these  conditions 
difficult . 

SEASONAL  RESULTS 

Both  models  described  in  this  report  are  functions  of  season.  In 
PROPHET,  the  seasonal  dependence  is  contained  in  the  LUF  and  MUF  predictions. 
For  HFBC84,  seasonal  dependence  is  found  in  the  MUF  and  in  the  ionospheric 
absorption  loss  term.  There  is  also  seasonal  dependence  in  the  auroral 
absorption  table,  as  discussed  above. 

To  investigate  the  accuracy  of  the  modeled  dependency  on  season  in  both 
models,  we  partitioned  the  data  bases  with  respect  to  the  four  seasons  and 
compared  predictions  with  data  for  each  case.  The  results  for  frequencies 
less  than  the  MUF  are  shown  in  figures  31  to  36. 
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Figure  23.  Average  absolute  relative  residual  as  a  function  of  circuit 
length  for  frequencies  less  than  or  equal  to  the  MUF. 
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Figure  24.  Correlation  coefficient  as  a  function  of  circuit  length  for 
frequencies  less  than  or  equal  to  the  MUF. 
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Figure  26.  Root  mean  square  residual  as  a  function  of  circuit  length  for 
frequencies  greater  than  the  MUF. 
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Figure  29.  Average  absolute  relative  residual  as  a  function  of  circuit 
length  for  frequencies  greater  than  the  MUF. 
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Figure  30.  Correlation  coefficient  as  a  function  of  circuit  length  for 
frequencies  greater  than  the  MUF. 
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In  figure  31  we  show  the  average  residual  for  both  models.  Each  model 
shows  a  peak  In  the  winter,  although  of  differing  sign.  Smaller  peaks  also 
occur  In  summer. 

On  the  other  hand,  figure  32  shows  that  the  dispersion  of  the  error,  rms 
residual,  is  much  larger  in  PROPHET  than  in  HFBC84,  with  a  peak  in  the  winter 
of  21  dB.  The  HFBC84  model  remains  approximately  constant  around  12  dB 
throughout  the  year. 

Similar  results  are  shown  in  figures  33  to  36  for  the  remaining 
statistical  parameters.  In  all  cases,  we  see  that  PROPHET  shows  large 
deviations  in  the  winter.  For  each  of  these  parameters,  HFBC84  performs  quite 
well  throughout  the  four  seasons,  with  relatively  high  correlation  as  shown  in 
figure  36. 

In  figures  37  to  42  we  show  a  companion  set  of  plots  for  frequencies 
exceeding  the  MUF.  Inspection  of  these  figures  shows,  in  each  case,  behavior 
similar  to  that  described  above.  The  PROPHET  model  shows  large  error  and 
variation  in  the  winter,  while  HFBC84  is  approximately  uniform  over  all 
seasons  in  most  cases. 

SUNSPOT  NUMBER  RESULTS 

In  HFBC84  the  ionospheric  absorption  loss  term,  equation  12,  shows 
dependence  on  12-month  running  mean  sunspot  number,  Dependence  on 

also  contained  in  the  long-path  model  through  the  absorption  equation  used  to 
determine  f  . 

Conversely,  the  PROPHET  model  contains  no  dependence  on  R^2’  Inspection 
of  equation  29  shows  the  QLOF  absorption  index  contains  no  sunspot-dependent 
component . 

To  determine  the  necessity  of  such  dependence  in  the  models,  we  performed 
data  comparison  based  on  the  parameter.  The  results  for  frequencies  less 
than  the  MUF  are  contained  in  figures  43  to  48. 

The  average  residual  shown  in  figure  43  indicates  there  is  no  strong 
dependence  on  R^2  bias  for  field  strength  predictions.  We  might  expect 

such  a  dependence  to  be  most  notable  for  large  values  of  R22’  however,  such 
behavior  is  not  indicated  in  this  figure. 
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for  frequencies  less  than  or  equal  to  the  MUF. 
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Figure  37.  Average  residual  as  a  function  of  season  for  frequencies 
greater  than  the  MUF. 
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Figure  38.  Root  mean  square  residual  as  a  function  of  season  for 
frequencies  greater  than  the  MUF. 
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Figure  39.  Average  relative  residual  as  a  function  of  season  for 
frequencies  greater  than  the  MUF. 
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Figure  40.  Root  mean  square  relative  residual  as  a  function  of  season 
for  frequencies  greater  than  the  MUF. 
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Figure  47.  Average  absolute  relative  residual  as  a  function  of  smoothed 
sunspot  number  for  frequencies  less  than  or  equal  to  the  MUF. 
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Figure  48.  Correlation  coefficient  as  a  function  of  smoothed  sunspot 
number  for  frequencies  less  than  or  equal  to  the  MUF. 
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Figure  44,  however,  shows  a  definite  dependence  on 
residual.  Performance  at  low  sunspot  numbers  is  very  poor  in  the  PROPHET 
model,  while  HFBC84  is  approximately  uniform  over  all  ranges. 

The  remaining  statistical  parameters  are  shown  in  figures  45  to  48.  These 
figures  show  that  the  strongest  dependence  for  field  strength  calculations  may 
be  for  low  sunspot  numbers  as,  in  each  case,  the  performance  of  PROPHET 
improves  for  high  sunspot  numbers.  In  each  case,  HFBC84  performs  quite  well 
with  only  minor  excursions  in  specific  sunspot  number  ranges. 

For  completeness,  we  present  in  figures  49  to  54  the  companion  plots  for 
frequencies  greater  than  the  MUF. 

MIDPATH  LOCAL  TIME  RESULTS 

Lastly,  to  determine  the  adequacy  of  the  time -dependent  aspects  of  each 
model,  we  looked  at  the  performance  as  a  function  of  midpath  local  time.  For 
long  paths  this  comparison  test  may  not  be  very  informative.  However,  for 
short  (1-hop)  paths  these  results  highlight  deficiencies  in  the  models. 

Figure  55  shows  the  average  bias  as  a  function  of  midpath  local  time  for 
frequencies  less  than  the  MUF.  The  figure  shows  a  diurnal  effect  in  PROPHET 
that  causes  overprediction  in  the  morning  hours  (-12  dB)  and  underprediction 
in  the  afternoon  and  night  (-12  dB) .  HFBC84  shows  quite  good  results 
throughout  the  day. 

Figures  56  to  59  show  no  clear  trends  in  the  other  statistical 
quantities.  Figure  57,  however,  shows  a  small  effect  similar  to  that  seen  in 
the  residual.  That  is,  a  tendency  to  overpredict  in  morning  hours  and  under¬ 
predict  in  afternoon  and  evening  hours. 

Figure  60  shows  the  correlation  coefficient.  In  this  figure  we  see  that 
PROPHET  shows  fair  performance  during  daylight  and  early  evening  hours. 
However,  its  nighttime  performance  is  again  quite  poor. 

Figures  61  to  66  show  a  companion  set  of  figures  for  frequencies  greater 
than  the  MUF.  In  this  set  of  figures,  neither  model  shows  any  clear  system¬ 
atic  trends.  For  the  most  part,  performance  is  equally  poor  in  both  cases. 

In  the  next  section  we  briefly  discuss  possible  explanations  for  the 
results  we  have  presented  here.  We  also  suggest  areas  of  research  to  improve 
these  results. 
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Figure  49.  Average  residual  as  a  function  of  smoothed  sunspot  number 
for  frequencies  greater  than  the  MUF. 
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Figure  50.  Root  mean  square  residual  as  a  function  of  smoothed 
sunspot  number  for  frequencies  greater  than  the  MUF. 


Figure  52.  Root  mean  square  relative  residual  as  a  function  of  smoothed 
sunspot  number  for  frequencies  greater  than  the  MUF. 
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Figure  53.  Average  absolute  relative  residual  as  a  function  of  smoothed 
sunspoi  number  for  frequencies  greater  than  the  MUF. 
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Figure  54,  Correlation  coefficient  as  a  function  of  smoothed  sunspot 
number  for  frequencies  greater  than  the  MUF. 
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Figure  57.  Average  relative  residual  as  a  function  of  midpath  local  time 
for  frequencies  less  than  or  equal  to  the  MUF. 
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Figure  58.  Root  mean  square  relative  residual  as  a  function  of  midpath 
local  time  for  frequencies  less  than  or  equal  to  the  MUF. 
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Figure  65.  Average  absolute  relative  residual  as  a  function  of  midpath 
local  time  for  frequencies  greater  than  the  MUF. 
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Figure  66.  Correlation  coefficient  as  a  function  of  midpath  local  time 
for  frequencies  greater  than  the  MUF, 
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DISCUSSION  AND  RECOMMENDATIONS 


In  the  previous  section  we  presented  the  results  of  a  comparison  of  the 
PROPHET  and  HFBC84  field  strength  predictions  to  measured  data.  In  this 
section  we  briefly  discuss  these  results  and  present  recommendations  for 
future  work  aimed  at  improving  them. 

First,  we  should  realize  that  comparison  of  predictions  of  monthly  median 
field  strengths  to  uncontrolled  observed  data,  as  we  have  done  here,  is,  in 
some  ways,  a  task  of  limited  validity.  Unlike  ionospheric  characteristics 
such  as  the  MUF  or  fQF2 ,  which  are  scaled  from  ionograms  in  the  same  way 
throughout  the  world,  procedures  used  to  obtain  field  strengths  vary  both  in 
method  and  in  quality. 

An  accurate  measurement  of  received  field  strength  requires  detailed 
knowledge  of  the  losses  inherent  in  both  receiving  and  transmitting  systems. 
This  includes  cable  and  antenna  feed  losses  at  both  ends  of  the  circuit.  In 
many  cases  such  detailed  knowledge  is  nonexistent. 

Furthermore,  to  apply  a  correct  antenna  gain  to  a  measured  field  strength 
we  must  know  which  mode  (or  modes)  we  are  receiving.  The  gains  are  then 
determined  by  prediction  programs  that  assign  take-off  angles  to  the  mode  ray. 
Any  errors  contained  in  the  antenna  gain  patterns  and  in  the  propagation 
prediction  programs  are  included  in  the  quoted  data  value.  These 
considerations  make  the  worldwide  collection  of  uncontrolled  data  for 
inclusion  in  a  data  base  a  questionable  process.  To  then  attempt  to  draw 
detailed  conclusions  from  comparison  of  predictions  to  sucl  data  is  a  risky 
proposition,  at  best. 

Aside  from  these  considerations,  the  overall  performance  of  the  models 
tested  in  this  report  is  quite  poor.  Specifically,  we  believe  that  the  over- 
the-MUF  prediction  models  require  further  investigations.  Uhile  we  usually 
work  below  the  MUF  in  any  operational  communication  system,  an  accurate 
prediction  of  usable  field  strength  at  frequencies  greater  than  the  MUF  could 
be  useful  in  certain  scenarios. 

Figures  1  to  6  show  that  both  models  perform  poorly  for  frequency-to-MUF 
ratios  greater  than  -1.5.  The  fact  that  the  average  relative  residual  and 
absolute  relative  residual  are  identical  in  this  range  shows  that  the  HF  ^f'84 
model  consistently  underpredicts  the  fiel  i  strength  in  these  cases.  This 
result  has  been  suspected  before  and  shows  up  clearly  here  (private 
communication,  David  B.  Sailors,  1986). 
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In  figure  67  we  show  a  comparison  of  the  current  CCIR  over-the-MUF  loss 

term,  L  ,  with  a  calculation  based  on  the  full  Wheeler  method,  equation  27. 
m 

Notice  that  the  CCIR  method  shows  a  bias  with  respect  to  the  full  Wheeler 
method.  We  recommend  a  calculation  based  on  this  method  be  implemented  in  the 
CCIR^  model  to  assess  its  accuracy  for  over-the-MUF  calculation. 


-  WHEELER  METHOD 

-  CCIR  METHOD  F/MUF 

Figure  67.  Comparison  of  full  Wheeler  method  and  current  CCIR 
method  for  over-the-MUF  loss. 

At  the  same  time,  the  PROPHET  model  similarly  shows  very  poor  results 
over  the  range  of  frequency-to-MUF  ratios  of  -1.5  and  higher.  This  implies 
that  the  determination  of  the  operational  MUF  using  the  HPF  factor,  which 
differs  from  the  FTZ  method,  is  not  adequate.  However,  the  poor  performance 
of  the  long-path  (FTZ)  model  in  HFBC84,  in  general,  shows  this  model  may  not 
be  capable  of  providing  accurate  global  predictions. 

Given  these  considerations,  we  recommend  that  future  work  on  the  PROPHET 
field  strength  model  be  geared  toward  the  Implementation  of  an  HFBC84-like 
model  for  paths  less  than  approximately  7000  km.  By  this  we  mean  a  prediction 
scheme  that  separately  models  the  various  loss  mechanisms  and  combines  them, 
as  in  equation  3,  to  produce  the  field  strength  nrediction. 
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A  desirable  feature  of  such  a  modeling  scheme  is  that  it  provides  a 
simple  method  for  addition  of  other  loss  mechanisms  that  may  be  developed. 
Terms  relating  to  these  additional  mechanisms  can  be  simply  added  to  equation 
3. 

In  the  same  way,  improved  models  for  the  existing  loss  mechanisms  can  be 
easily  inserted  and  tested.  For  example,  the  PROPHET  ionospheric  absorption 
term  used  in  QLOF  and  in  the  determination  of  f^  in  the  field  strength  model 
is  quite  accurate,  we  believe.  It  provides  excellent  performance  in  LUF 
predictions  and  may  be  preferable  to  the  expression  used  in  HFBC84. 

Based  on  the  results  of  the  comparison,  we  recommend  using  the  current 
PROPHET  implementation  of  the  FTZ  model  for  paths  greater  than  7000  km.  From 
figure  19  we  see  that,  except  for  the  very  long  paths  where  both  models  show 
poor  results,  PROPHET  does  quite  well  in  terms  of  the  average  bias  for  paths 
up  to  -10,000  km.  Similarly,  the  rms  residual,  figure  20,  for  the  PROPHET 
model  is  as  good  as,  or  better  than,  HFBC84  for  the  same  path  length  ranges. 
Further  work  is  needed  at  the  very  long  path  lengths  where  the  problem  is  much 
more  difficult  due  to  multimode  propagation  and  scattering  that  makes 
prediction  difficult. 

In  summary,  this  report  has  highlighted  the  good  and  bad  aspects  of  both 
prediction  models.  A  combination  of  the  best  aspects  of  both  into  one  model 
should  lead  to  a  prediction  scheme  which  is  fast,  accurate,  and  will  provide  a 
reliable  input  into  higher  level  models. 
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APPENDIX  A 


OVERALL  RESULTS 


This  appendix  presents  overall  results  in  tabular  form  for  the  cases 
investigated  in  this  report.  In  each  table,  the  average  residual,  rms 
residual,  average  relative  residual,  rms  relative  residual,  average  absolute 
relative  residual,  and  correlation  are  given  for  each  condition.  The  values 
given  in  the  table  are  in  the  form  HFB84/PR0PHET . 


TABLES 


A-1.  Summary  of  results  as  a  function  of  a  frequency- to -MUF  ratio  for 

entire  data  base  .  A-2 

A- 2.  Summary  of  results  as  a  function  of  circuit  for  entire  data  base  A- 3 

A- 3.  Summary  of  results  as  a  function  of  circuit  length  for  entire 

data  base .  A- 5 

A-4.  Summary  of  results  as  a  function  of  season  for  entire  data  base  A-6 

A- 5.  Summary  of  results  as  a  function  of  smoothed  sunspot  number  for 

entire  data  base  .  A-6 

A-6.  Summary  of  results  as  a  function  of  midpath  local  time  for 

entire  data  base  .  A-7 
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Table  A-1.  Siiitimary  of  results  as  a  function  of  frequency-to-MUF  ratio  for 
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Table  A-4,  Summary  of  results  as  a  function  of  season  for  entire  data  base. 
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